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The programming new e + e~ collider with high luminosity shall provide another useful platform 
to study the properties of the doubly heavy B c meson in addition to the hadronic colliders as 
LHC and TEVATRON. Under the 'New Trace Amplitude Approach', we calculate the production 
of the spin-singlet B c and the spin-triplet B* mesons through the Z° boson decays, where 
uncertainties for the production are also discussed. Our results show IVs.) = 81. 4+4Q 2 j 1 KeV 
and r (3Sl) = 116.4±^ 3 g 9 KeV, where the errors are caused by varying mt and m c within their 
reasonable regions. 
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The B c meson is a double heavy quark-antiquark 
bound state and carries flavors explicitly. Since its first 
discovery at TEVATRON 0, B c physics is attracting 
more and wide interests. Recently, many progresses have 
been made for the hadronic production of B c meson at 
high energy colliders as LHC and TEVATRON. A com- 
puter program BCVEGPY for the direct hadronicpro- 
duction of B c meson has been presented in Rcfs.JJQ- 
And it has been found that the indirect production of B c 
via top quark decays can also provide useful information 
on B c meson 0-0] • 

Comparing with the hadronic colliders, an e + e~ col- 
lider has its own advantages, mainly because of its lower 
background. As for the previous LEP-I experiment, no 
B c events have been found due to its lower collision en- 
ergy and low luminosity However, if the luminosity 
of the e + e~ collider can be raised up to C oc 10 34 cm _2 .s _1 
or even higher as programmed by the Internal Linear Col- 
lider (ILC) H3, then there might have enough events. 
Moreover, if the e + e~ collider further runs at the Z°- 
boson energy, the resonance effects at the Z° peak may 
raise the production rate up to several orders. It has 
been estimated by Ref.[ll[ that more than 10 9 ~ 10 Z°- 
events can be produced at ILC per year, which is about 
3^4 orders higher than that collected by LEP-I. Such 
a high luminosity collider is called as GigaZ [ll| or a 
Z-factory [13]. Then it will open new opportunities not 
only for high precision physics in the electro- weak sector, 
but also for the hadron physics. 

The production of B c through Z° decays has been 
studied in Refs [1, [t| [l3| with various methods. Since 
the process is very complicated, it would be helpful to 
have a cross check of these results. Furthermore, con- 
sidering the forthcoming Z-factory, it may be interesting 
to know the theoretical uncertainties in estimating of B c 
production. 

For the purpose, we need to calculate the process 
Z° — s- Bc*^ +6+c, whose Feynman diagrams are shown in 
Fig- (H|). According to the NRQCD factorization formula 
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FIG. 1: Feynman diagrams for the process Z°(k) — > 
Bi*\ q3 ) + b(q 2 ) + c(qi). 

14], the decay width for the process Z° —t Be +b + c 
can be written in the following factorization form: 

dT = ^2df(Z° -> cb[n} + b + c)(0 H ), (1) 

n 

where the matrix element (O h (n)) is proportional to the 
inclusive transition probability of the pcrturbative state 
cb[n] into the bound states of B c . As for the two color- 
singlet S-wave states c6[ 1 5o] and c6[ 3 Si], their matrix 
elements (O h (n)) are related with the Bethe-Salpeter 
wave function at the ori gin that can be determined by 
the potential model [l5l-|20|. dt(Z° -> cb[n] + b + c) 
stands for the short-distance decay width, i.e. 

dT(Z° cb[n] + b + c) = J-^|M| 2 d$3, (2) 

where ^ means we need to average over the spin states 
of initial particles and to sum over the color and spin of 
all the final particles. And in Z° rest frame, the three- 
particle phase space can be written as 
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The hard scattering amplitude for the process Z°(k) —> 
Bc*\q 3 ) + b(q2) + c(q±) can be written as: 

4 

iM = Cu 8 (q 2 ) }^ A n v s >(gx), (3) 



n=l 



where C= sin g m C ose m x 3~7r Tne g amma structure A„ 
(n = 1, • • 4) corresponds to the four Feynman diagrams 
in Fig.([T]). More explicitly, A n can be written as 
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and q is the relative momentum between the two con- 
stitute quarks of cfr-quarkonium. In the nonrelativistic 
approximation, the S'-wave projector xf 3 (<z) takes the 
following form 



x£f"(?) 



Am b m c 



($ 31 -™>b) (ary 5 it 32 + m c) 

(8) 

where a = 1(0) and j3 = 0(1) for S = 0(1) meson respec- 
tively. s s (q 3 ) is the polarization vector, 931 and 932 are 
the momenta of the two constitute quarks of the meson, 



931 



q 3 + q and q 32 



-q 3 - q, 



(9) 



where mj c = m b + m c is implicitly adopted to ensure the 
gauge invariance of the hard scattering amplitude. 

By using the conventional trace technique, we need 
to derive the squared amplitude, which is very compli- 
cated and lengthy for the present case. To derive ana- 
lytical expression for the process Z° — > Be + b + c and 
to make its form simpler as much as possible, we adopt 
the 'new trace amplitude approach' suggested and devel- 
oped by Refs. 0, H[ to do our calculation. Under the ap- 
proach, we first arrange each of the four amplitudes listed 
in Eqs. (|4l5l6l7[) into four orthogonal sub-amplitudes ac- 
cording to the four spin combinations of the outgoing 
6-quark and c-antiquark, and then do the trace of the 
Dirac-7 matrix strings at the amplitude level by prop- 
erly dealing with the massive spinors, which results in 
explicit series over several independent Lorentz-invariant 
structures. And then, our task left is to determine the co- 
efficients of these Lorentz-invariant structures. To make 



the paper more compact, we present the detailed formu- 
lae for dealing with the process in Appendices A and B, 
where Appendix A gives the phase-space integration for- 
mulae and Appendix B gives the 'new trace amplitude 
approach', which presents all the necessary coefficients 
for the Lorentz-invariant structures. 

As a cross check of the present obtained results, by 
taking the same parameters, we can obtain consistent 
numerical results as that of Ref . [§| within reasonable nu- 
merical errors 1 . 

In doing the numerical calculation, we take mz = 
91.1876 GeV and a s (m z ) = 0.1176 [p. To be consistent 
with the present leading-order calculation, we adopt the 
leading-order a s running, and by taking the normaliza- 
tion scale to be 2m c , which leads to a s (2m c ) — 0.212. 
The two constitute quark masses are taken as m b — 4.90 
GeV and m c — 1.50 GeV. With the above parame- 
ter values, it can be found that the total decay width 
r ( i So) = 81.4 KeV and r ( 3 Sl) = 116.4 KeV. 

The differential distributions of the invariant masses 
si and s 2 , i.e. dT/dsi and dT/ds2 are shown in Fig. ([2]), 
where Si = (qi+q 3 ) 2 and s 2 = {qi + q 2 ) 2 ■ And the differ- 
ential distributions of cos 6*13 and cos $23, i.e. dT /d cos 8i 3 
and dT/dcos8 23 are shown in Fig.®, where 8i 3 is the 
angle between qi and q 3 , and #23 is the angle between 
q 2 and q 3 respectively. It can be found that the largest 
differential decay width of dT/d cos 6*13 is achieved when 
#13 = 0°, i.e. the (c6)-quarkonium and c-quark moving in 
the same direction. While the largest differential decay 
width of dT/d cos 6*23 is achieved when #23 = 180°, i.e. 
the (c6)-quarkonium and 6-quark moving back to back. 

Next, it would be interesting to show the theoreti- 
cal uncertainties for the production. Main uncertainty 
sources include the matrix elements (or the wavefunc- 
tion at the origin of the binding system („)(0)|), the 
renormalization scale fiR, the constitute quark masses 
m b and m c . \tp R M (0)| and «(//#) are overall parame- 
ters for the present case, and their uncertainties can be 
easily figured out. For example, one can set (1r to be 
2m c or 2m b , since the intermediate gluon as shown in 
Fig.(fT]) should be hard enough so as to produce a cc- 
quark pair or a 66-quark pair, which inversely ensures 
the pQCD applicability of the process. By setting these 
two scales to calculate the process, we obtain the ratio 
r MJJ =2m i) /r A1R= 2m c oc a 2 (2m b ) / a 2 (2m c ) ~ 0.67. In the 
following discussion, we fix (1r — 2m c and |^ B («)(0)| = 

0.361 GeV 3 / 2 go|. 

For clarity, we present the uncertainties of m c and m b 
in 'a factorizable way'. When focussing on the uncertain- 
ties from m c , we let it be a basic 'input' parameter vary- 
ing in a possible range m c = 1.50 ±0.30 GeV with all the 
other factors, including the 6-quark mass and etc. being 
fixed to their center values. Similarly, when discussing 



1 There are some typos in the formulae listed in the Appendix of 
Ref. the right ones are presented in the present Appendix B. 
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FIG. 3: Differential decay width dT/d cos 613 (Left) and dr/dcos#23 (Right) for the process Z — > B * +b + c, where the solid 
and the dashed lines are for B c and £>*, respectively. 



m c (GeV) 


1.20 


1.50 


1.80 


r ( i Sn) (KeV) 


183.5 


81.4 


42.2 


r (3Sl) (KeV) 


280.1 


116.4 


57.1 



TABLE I: Decay width for the production of Be through Z' 
decay with varying m c , where mt is fixed to be 4.9 GeV. 



rn b ( GeV) 


4.50 


4.90 


5.30 


r ( i Sn) (KeV) 


82.1 


81.4 


71.0 


r ( 3 Sl) (KeV) 


114.1 


116.4 


95.6 



TABLE II: Decay width for the production of through 
Z° decay with varying mj, where m c is fixed to be 1.5 GeV. 



the uncertainty caused by nib, we vary the 6-quark mass 
nib within the region of nib — 4.90 ± 0.40 GeV. 

The decay width for the production of Be through Z° 
decay with varying m c or nib are presented in TAB U and 
TAB [II] It shows that the decay width is more sensitive 
to m c , which decreases with the increment of m c . 



By adding these two uncertainties caused by nib and 
m c in quadrature, we obtain 

I>s ) = 81.*?™ KcV (10) 

and 

r ( 3 Sl) = 116.4^ 3 8 9 KeV. (11) 

The shaded bands in Figs. (14l5[) show the corresponding 
uncertainty more clearly, where the contributions from 
1 S'o and 3 Si are summed up. The center solid line is 
for m c — 1.5GeV and nib = 4.9GeV, the upper edge of 
the band is obtained by setting m c = 1.2GeV and nib — 
5.3GeV, while the lower edge of the band is obtained by 
setting m c — 1.8GeV and nib = 4.5GeV. 

As a summary: by using the 'new trace amplitude ap- 
proach', we calculate the B c production process, Z° —> 
B c ' + b + c. The decay widths together with their uncer- 
tainties caused by the b and c quark masses are rVi,s ) = 
SlAt^ 1 KeV and r (3gl) = 116.4^| 3 g 9 KeV, where the 
errors are caused by varying nib and m c within their rea- 
sonable regions m c g [1.2,1.8] GeV and m c G [4.5,5.3] 
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FIG. 4: Uncertainties of differential decay width dY/dsi (Left) and dY /ds2 (Right) for Z° — s> B c + b + c, where the contributions 
from and 3 Si are summed up. 
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FIG. 5: Uncertainties of differential decay width dY /d cos #13 (Left) and (iT/a! cos #23 (Right) for Z° — > B c + b + c, where the 
contributions from 1 So and 3 .Si are summed up. 



GeV. Further more, the differential decay width for sx t 2 
and cos #13^3 together their uncertainties are drawn in 
Figs. (|4l5|) . Considering the advantage of the clean envi- 
ronment in e + e~ collider, it will provide another useful 
platform in studying the B c production. 
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Appendix A: Formulae for the phase space 
integration 

The decay width of the process Z°(k) — >• Bl*\q 3 ) + 
Hq?) + c(gi) is proportional to the following phase space, 



where k = (fc°,fc) = (fc°, fc 1 , fc 2 , fc 3 ), q f = (q° f ,q f ) = 
(q®, q l j,q 2 j,q^). Furthermore, in the rest frame of Z° bo- 
son (fc° — mz), we have 
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« TTTTTTo 7 Ai A 2 A 3 %i - to 2 )<5(<z 2 2 - m 2 )<K? 3 2 - m| c )%°)% 2 )% 3 ) 



dsids 2 (27r) 5 (2m z ) 

x<5(si - (gi + g 3 ) 2 )<5(s2 - (<?i + q 2 ) 2 )5{k - q x - q 2 - (73) 

OC ¥ ^-d l Q i d i q 3 5{{k -q 2 - q 3 ) 2 - m 2 c )5(q 2 - m 2 )S(q 2 - m 2 Br )9(k° - q° 2 - q° 3 ) 

xO(q° 2 )6(q 3 )5( Sl - (k - q 2 ) 2 )S(s 2 - (k - q 3 ) 2 ) 
OC ¥ ^ 3 d 3 q 2 d 3 q 3 6(qf - f 2 - m 2 c )S(qf ~ <fs ~ m 2 B J0(k" q° q") 

xO(q° 2 )0(q 3 )S( Sl + m% c - m 2 c - 2m z q° 3 + 2q° 2 q° 3 - 2q 2 ■ q 3 ) 
oc Ml^-dfia sin 9 23 d9 23 dcb 23 9(k - q° 2 - q^9(q" 2 )9(q° 3 ) 

xS(s 1 +s 2 -m%- m 2 c + 2qlql - 2\q 2 \ ■ |<7 3 | cos <9 23 ) 
« 2%t| 0(fc ° ~ q2 ^)0^)0{ql)9{X), (A2) 



where q 2 = m z+ m b - s i an( j _ m z+™B c s 2 ^ ^ _ With all these step function above, we can get the intc- 



2m z *3 2m z 

q 2 2 —m 2 and \q 3 \ — \J q^ — m 2 c . The step function 
9{X) is determined by ensuring | cos# 23 | < 1, where 

St + s 2 - m| - to 2 + 2^9° 

COS«23 = ^TZTTTZri • 



gration ranges: 



(s 2 - to| + m|J (s 2 - to 2 + to 2 ) + J77 (s 2 , to|, m|J 77 (s 2 , to 2 , to 2 ) 

s r n = toc+^l 2^ (A3) 

(s 2 - m% + m% ) (s 2 - to 2 + to 2 ) - Jjj (s 2 , m|, m^J 7/ (s 2 , m 2 , to 2 ) 

^ ax - to' + "»| e — I (A4) 

s f n = (m c + m b ) 2 (A5) 

= (toz-tobJ 2 (A6) 



where 7/(2;, y,z) = (x — y — z) 2 — Ayz. 

Furthermore, we can obtain the cos 6* 2 3 distribution 



and there are two s 2 in different range of cos #23 and S\, 
i.e. 



dT 
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dsidcos9 23 2 7 7t 3 to| 
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where the extra Jacobian 



J = 



1*32 1 \fo\ 



I _ ji_ , l92|(m|+m| e -s 2 ) 

™z m z ^ m 4 Bc -2(;ml + s 2 )ml a + (rnl-s 2 r 



cos 9 23 
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|92| 2 ( to b c + m z 2 ) cos 2 6*23 - (<?2 ~ ™>z)[m z {si + q 2 m z ) + q 2 m 2 Bc - m \- m z" 



\q 2 \ 2 cos 2 6»23 - (<?2 ~ TO ^) 2 

±mz|g5| cosfl23[mi B - 2m 2 3c (m 2 + 2(m z - q a 2 ) 2 - si - 2|<f 2 | 2 cos 2 23 ) + (m 2 . - Sl ) 2 }? L (A9) 



where s 2 is obtained when cos6*23 G [0,-1] and si G or cos#23 G [0,-1] and si € [si m i n [cos #23], si max]- The 
[simin[cos023],simin[cos023 = 0]]. And is obtained 0(X) function determines the boundary of sj.: 
when cos 23 G [1,0] and si G [si min [cos 6» 23 = 0],si max ] 



simax = (m z -m b ) 2 (A10) 



„ , m^(cos 2 6» 23 - l)m| + m|(m 2 + m| cos 2 6» 23 ) + m Bc mzVY /A11 . 

SlminCOS023 = ? Tfl TV~^2 \ 2 ( A11 ) 

(cos 2 023 - l)m B + m z z 
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(( m Bc 
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The distribution for cos #13 can be obtained in a similar 
way. 



Appendix B: Amplitude of the process 

Z°(k)^ Bi* ) (q3) + b(q 2 ) + c(q 1 ) 

The amplitude M of the process Z°(k) -> Bi*\q 3 ) 
6(32) + c(gi) has the general structure 



M = u s (g 2 )^w s '(9i) 5 



(Bl) 



where A can be read from Eqs.Q-©. 

To derive analytical expression for the process and to 
make its form simpler as much as possible, we adopt the 
'new trace amplitude approach' suggested by Refs.0, [|| 
to do our calculation. Detailed process of the approach 
can be found in Refs.fl, [||, and here, we shall only list 
our main results. 

After summing up the spin states, the square of the 
amplitude can be divided into four parts, 



\M\ 



\Mx\ 



IM2I 2 



Wh\ 2 



|M 4 | 2 , 



(B2) 



where by introducing a light-like momentum kg and a 
spacelike vector k\ that satisfies the relations, k\-k\ = — 1 
and kg ■ k\ = 0, the four amplitude Mi can be written as 



Mi 
M 2 
M 3 



N 

71 

N 

V2 
N 

71 



Tr 
Tr 
Tr 



{fa - "lc)^o(^ 2 + mb ) A 

{fa ~ m c )7 5 ^ (^ 2 + m b )A 
{fa ~ m c )loW^ 2 + m b )A 



and 



N 

M 4 = ^Tr 
y/2 



{fa - "lc)75^l^o(^ 2 + m b) A 



(B3) 



where N — l/-y/4(/c ■ qi){k ■ q 2 ) is the normalization 
constant. k and ki are arbitrary momenta, and in order 
to write down M n as explicitly and simply as possible: 
1) We set k = q 2 — aqi, where the coefficient a is deter- 
mined by the requirement that fco be a lightlike vector: 



qi ■ 92 ± \/{qi ■ <72) 2 



— 5 — 9 

mimi 



(B4) 
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2) We set fcf = iN a e^ up<y qi v k p q 2ai where AT ensures k\ ■ 
k\ = — 1. It is found that ^ can be expressed as, 

^ 1 = JVo75 91 • fc^ 2 + j^fc • g 2 - 9i • q2$ - ify 2 ■ (B5) 
And then the resultant M; can be simplified as: 

M 1 = L x xTr[(^-m c )(^ 2 +m 6 )A], (B6) 

M 2 = L 2 x Tr[(^ - m c )7 5 (j£ 2 + m b )A], (B7) 

M 3 = M 3 > - N [m b (qi ■ k) + m c (q 2 ■ k)]M 2 , (B8) 

M 4 = M 4 < + N [m b ( qi ■ k) - m c (q 2 ■ k)]M u (B9) 



and 



where 



M 3 , = ^Tr 
4L 2 



Ma, 



4Li 



Tr 



!i - m c ) 75 % 2 + m 6 )^lj , (BfO) 
- m c )#(^ 2 + m b )A] . (Bff) 



Furthermore, the amplitudes Mj can be expanded over 
some basic Lorentz structures: 

m 

M i (n)=Y, A i(n)B j {n){i = l-A) (B12) 
i=i 

and 

m 

Mji(n) = ^]4'(n)S j (n) (i' = 3,4) (B13) 
i=i 

where m is the number of basic Lorentz structure Bj (n) , 
whose value dependents on the (c6)-quarkonium_state n: 
e.g. m = 3 for n = (c&)[ 1 5o]i, m — 12 for n — (cb)[ 3 Si]i. 
As for A?(n) and Alj(n), they can be expressed by 



x = 93 • k/m% = 1^2 ( m B c + m z - s 2), 
V = 92 • k/m% = ^r{ m b + m%- si), 
z = q 1 -k/m% = ^r( m z +m 2 c - s 3 ), 
w = 93 • 92/m| = T^T ( s 3 - m| c - mg), 



v = 93 • 9iMl : 
w = 9i • W TO I 



2m| 



2m' 



(si - m| - m 2 ), 



-(s 2 -m 6 -m c ), 



where si = (91+93) , « 2 = (91 + 92) , and s 3 = (92 + 93) , 
which satisfy the relation: si + S2 + S3 = m| + to 2 + m 2 + 
m 2 ^. And the short notations for the denominators are 



di 
d 2 
d 3 

G?4 



{q 2 -kf-ml (931 +9l) 2 ' 

1 1 
(k - 932) 2 - ml (931 +91) 2 ' 

1 1 
(932 + 92) 2 (931 - k) 2 - ml ' 

1 1 
(932 + 92) 2 (93 + 92) 2 - ml' 



Furthermore, the following relations are useful to short 
the expressions: 

u + v + r 2 = x, w + u + r\ = y, w + v + r\ = z. 



A](n) = Af(n) - N [m b ( qi ■ k) + m c (q 2 • k)]A 2 {n), 
A$(n) = Af(n) + N [m b ( qi ■ k) - m c (q 2 ■ k)]A)(n). 

The explicit expression for A^' 2 (n) and ' 4 (n) of each 
state shall be listed in the following subsections. 

To shorten the notation, we set T b = j — isin 2 ^,,, and 



j — |sin 2 0„,. And, wc define some dimensionless 
parameters 



n = 



m b 
mz 1 



r 2 = 



m c 
m z ' 



r 3 = 



m Br 
m z 



1. Coefficients for the production of B, 

;re are 3 basic 
B c ( 1 5q), which are 



There are 3 basic Lorentz structures Bj for the case of 



93-e(fc) D 92-e(fc) „ i n 
Bi = , B 2 = , B 3 = —e(k, 93,92, e(k)), 



m z 



m z 



mi 



where e(k,q 3 ,q 2 ,e(k)) = s> il,p ' J k p q 3lJ q 2p e rT (k). The values 
of the coefheients A 1 - and A? are 



A\ = Lim ^ - 2rxr 3 ) - (2n - r 3 )y)di + (n(r 3 + 2r 3 u - 2r lX ) - rly)d 2 + (n(2r 2 (a; - u) + r 3 (y - 1)) 

+r 2 (r 3 y - 2r 2 u))d 3 + {r\ + 2(r 2 - 2r 3 )r 2 + (r 2 - r\ - y)n + (r 2 - 2r 3 )(2w - y))d A ) , (B14) 



A\ = -^—=—((r 3 (x-2u) + r 3 (Au-2x-Ay + 2))d 1 + (r 3 x-2r 1 r 2 )d 2 



+ {2r 2 r 2 - r 3 x)d 3 + {-2 ri r 2 + 2r 2 r 2 - 2ur 3 - xr 3 )d 4 ), (B15) 
A\ = -AL l m z 7/2 ^T b {d l +d 2 + d 3 + d 4 ), (B16) 

A\ = ° (T b (-2y 2 + y- nr 3 )r 3 rf 1 + T b {-2xr\ + (r 3 (4x + Ay - 3) - 2r 2 a;)r 2 + (r 2 r 3 -2u + 2(r 2 - 

i 2 r 3 7 

2r 2 r 3 + 2u)y)n + r 3 (2u + {-2x - 2y + l)y))d 2 - T c ((r 3 - 2r 2 x)r 2 + r 2 (r 3 (4x + 4y - 3) - 2r 2 a;)n 

-2r 3 y 2 - 2r 2 w + (-4r 3 r 2 + 2r 2 r 2 + 4w 2 + r 3 )y)d 3 + T c (nr 3 (l - 2x) + (2r 3 (r 3 - 3n) - 2u)y)r 3 d 4 ), (B17) 

A 2 ' = f " (Tb^xrl + 2(-2xr 3 - 2yr 3 + r 3 + r 2 x)r 4 + (-r 2 + 2r 2 r 3 - 2u + x) x (2y - l))di 
L 2 y/r 3 

-T b (r\ + x(-2x -2y + l))d 2 - T c (r 2 3 + x(2y - \))d 3 + T c (xr 2 + 2(-2xr 3 - 2yr 3 + r 3 + r 2 x)r x 

+r 2 2 x + 2r 2 r 3 (2y - 1) - (r 2 + 2u)(x + 2y - l))d 4 ), (B18) 

A 3 3 ' = mz ' ^ {{2nr 3 + 2y- l)di + (2nr 3 - 2x - 2y + l)d 2 - (2r 2 r 3 + 2y - l)d 3 + (2r 2 - 4r 2 r 3 + 2u)d 4 ) (B19) 
AL 2 y/r 3 

The values of the coefficients A 2 and A* are 

A\ = - 4L2 ^L /2 (T b ((2n + r 3 )y - n)di + T 6 ((5n - r 2 )r 3 y - r^nr 2 + (4n(n - r 2 ) + 4w 

+l)r 3 + 2riu - 2r 2 w - 2rix))d 2 + (rir 2 r 3 T c (6ri - 2r 2 ) + r 3 T c ri + (4r 2 ri + 2r 2 r 3 )T c u 
-2r 2 T c xri - r 3 T c yr x - 2>r 2 r 3 T c y)d 3 + (-r 3 T c r\ + Ar\T c r\ + r 3 T c n - Ar 2 r\T c r x 

+r 2 r 3 T cri + r 3 T c y ri + Ar 2 T c u - 2r 2 r 3 T c u - 2r 2 T c y + r 2 r 3 T c y)d 4 ), (B20) 

AT ^ 7/2 

A 2 = — r = — (T 6 ((6n - 2r 2 )rir 3 + (4u - 2x - Ay + 2)r 3 + 2r x u - 2r 2 u - 3r lX + r 2 x)d 1 

+ (r 2 - n)T b xd 2 + (n - r 2 )T c xd 3 - T c (-2r 3 ((n - r 2 ) 2 - r 2 r 3 ) - 2(n - r 2 )u - r 3 x)d 4 ), (B21) 
^ = L2m !. 7/2 (r 3 di + (r 2 - n)da + (r 2 - n)d 3 - r 3 d 4 ), (B22) 

- ^^-^((r 1 r 3 + y)(2y-l)rf 1 + (2a;r3-(r 3 + 2r 2 x)r 2 + (r 2 r 3 + 2 M -4 U y)r 1 
-r 3 (2u + (-2x -2y + l)y))d 2 + ((r 3 - 2r 2 x)r 2 + r 2 (2r 2 x - r 3 )n 

+ (2y - l)(2r 2 u - r 3 y))d 3 + (nr 3 (2x - 1) - (-4nr 3 - 2u)j/)d 4 ), (B23) 

A% = m A Z T ' ^ ((r 3 (2n - 2r 2 + r 3 ) + 2u - x)(2y - l)di + (r 2 + 2n(.x + 2y - l)r 3 + .t 
4i>iV r 3 

-2a;(a; + y))d 2 + (r 2 - 2r 2 (a; + 2y - l)r 3 + x(2y - l))d 3 + (-xr 2 + 2r 3 (2x + 2y- l)n 

+r 2 a; + (r 2 + 2?i)(.x + 2y - 1) - 2r 2 r 3 (2a; + 2y - l))d 4 ), (B24) 

A A 3 = _ N ° mz ^_ 2 ( Tb (2y - l)d! + T b (-2x - 2y + l)d 2 + (T e - 2T c y)d 3 + T e (r 2 - r 2 + r 2 + 2 M )d 4 ) (B25) 



I 

Coefficients for B c * R fc ' e (g3)g3 ; e{k) k ■ e(q 3 )q 2 ■ e(k) 

t>5 = 2 ' -°6 = 2 : 



'z 



There are 12 basic Lorentz structures Bj for the case R q 2 • e(93)?3 • R _ <Z2 • c(k)q 2 ■ e(3 3 ) 

of B* ( 3 ^), which are ^ 7 _ ^ ' ^ 8 ~ ' 



1 

d /,s / n d * n /; \ / \\ B 9 = —e(k,q 3 ,q 2 ,e(k))k-e(q 3 ), 

B i = ^2 = — £(fc,<?3,e(fc),e((? 3 )), m| 



m| 



d * n /7 \ / \\ r? * / n \ ( \\ Bl o = — r £ ( fc >93,<72,e(<?3))<73 • 

S 3 = — e(k,q 2 ,e(k),e(q 3 )), B 4 = —e(q 3 ,q 2 ,e{k),e(q 3 )), m% 

Tfl y 777 7 
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fin = -\e(k, q 3 , q 2 , e(q 3 ))q 2 ■ e(fc), The values of the coefficients A] and Af are 

B12 = —e(k,q 3 ,q 2l e(k))q 2 -e(q 3 ). 
m z 



A\ = 4Ll ^ /2 (Tb{r 2 y + n(x + y- l))r 3 di + T t (-2ay? + (r 3 (a; - 1) - 2r 2 x)n 

~2ux + r\y + 2xy)d 2 - T c (2xr\ + (2r 2 x + r 3 (x - l))n + 2ux + r§j/ - 2xy)d 3 

-T c {r\ -(r 2 2 -r 2 -2u + x + y)n + r 2 y)r 3 d 4 ), (B26) 

T 7/2 

A\ = im f_ {-nr 3 di + (n(-2n - r 3 ) - 4u + 4 ? y)d 2 + (2r? + r 3 n + Au - Ay)d 3 - nr 3 d 4 ), (B27) 

V^3 



r ^ 7/2 

A\ = 1 £. ((n - r 2 )r 3 d 4 + (3r§ - 2x)d 2 - (r§ - 2x)rf 3 - (n - r 2 )r 3 rf 4 ), (B28) 

22 7im ^7/2 

^4 = -== — (r"ir- 3 rfi + (nr 3 + x — l)d 2 + (r 2 r 3 - x + l)d 3 + r 2 r 3 d 4 ), (B29) 



A\ = 4Ll ™| ' (-nr 3 T b di + T b (2r 2 + ran + 2m - 2y)d 2 + T c {2r 2 + r 3 n + 2u - 2y)d 3 - nr 3 T c d 4 ), (B30) 



A\ = -ALxmz 7 ' 2 ^ + r 2 )V^T b di + r 3 3 /2 T b d 2 - r 3 3 /2 T c d 3 + (n - r 2 )^T c d 4 ), (B31) 
A\ = 8L imz 7 ^(^r lv ^T b d 2 +r 2 ^T c d 3 + rl /2 T c d 4 ), (B32) 
A\ = 8L 1 m z 7 ' 2 rl /2 {T b d 1 +T c d 4 ) 1 (B33) 

Al = A\^ 2L ^'\ d 2 -d 3 \ (B34) 

A\ x = A\ 2 = 0, (B35) 

AX = - m . Z T ' '^ (r 3 ((2x + 2y - l)r\ + (r 2 - 2r 2 y)n +u + y- 2y(x + y))di + (-2xrf + (r 3 (Ax + Ay - 3) 
AL 2y /r 3 

-2r 2 x)r\ + (r 2 r 3 - 2(u + (x - l)x) + 2{r 2 3 - 2r 2 r 3 + 2u - x)y)r 1 + 2r 2 xy + r 3 (-2y 2 - 2xy + y + u))d 2 
-((r 3 - 2r 2 x)r^ + {-2xr\ + r 3 (Ax + Ay - 3)r 2 - 2x(x + y- l))n - 2r 2 u - Ar\r 3 y + 2r 2 (r\ + 2u + x)y 
+r 3 {-2y 2 - 2xy + y + u))d 3 + r 3 (-u - n(n + r 2 (2x - 1)) + ((n - r 2 f + r\ + 2u)y)d 4 ), (B36) 

At = 2m f' 2 ^ (T b (n(x - 1) + (n - r 2 )y)d 2 - r^yd^ + T c ( n (x - 1) + (n - r 2 )y)d 3 - r 3 T c (r? + u)d 4 ), (B37) 
At = ^ T N °^ (T b (l - 2y)d x + T b (2x + 2y - l)d 2 + T c (2x + 2y - l)d 3 + 2r 2 r 3 T c d 4 ), (B38) 
At - - mz9/ r Noy ^ (Thidi + d 2 ) + T c d 3 + T c (2x + 2y- l)d 4 ), (B39) 

-^2 

= ^| ^ (r 3 (r^ - y)d 1 + (r 2 y -n(u + x + y- l))d 2 + (r 3 y - r 2 (u + y) + r 1 (x + y- l))d 3 + r 3 {r\ + v£)BA))) 

^^2V r 3 

At = mZ / A ^°^ (-(2u + 2y- l)d 4 + (2nr 3 -2x-2y + l)d 2 + (2r 2 r 3 + 2y - l)d 3 + 2r 2 r 3 d 4 ), (B41) 
AL 2 

A 3 7 ' = - mz '^ (-r 3 rfi + (r 3 - 2nx)d 2 + (r 3 - 2r 2 a;)rf 3 - r 3 (2z + Ay - l)d 4 ), (B42) 
AL 2y /r 3 

4 = mz9/2 " 0Vr3 (x + 2y-l)(d 1+di ), (B43) 

ZL 2 

4 = - 2 " lz ' / ; Ar °^ r^ 3 , ,B44) 
L 2 



A$ = (r 3 (2n(r 2 - n) + 2y - l)rfi + (-2nrg - 2yr 3 + r 3 + 2(n - r 2 )x)d : 
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< = A?! = 2m T z9/ ^ ( ri T b d 2 r 2 T c d 3 ), (B45) 

^2 

The values of the coefficients A 2 - and A A - are 

A\ = 2 1. ( r 3(2r? - 2r 2 r 2 + (2u - x - 3y + l)n + r 2 y)rfi + (r 3 n(6r 2 - 2nr 2 ) 

+2(u — 2a;)r 2 + (2r 2 (a; — m) + r 3 (4u — a; — 5y + l))ri — 2wa; + r 2 r 3 y + 2xy)d 2 
-((4r 2 r 3 - 2x)r\ + (2r 2 r\ - (Ar\ + x + y - l)r 3 + 2r 2 u)r 1 - 2r\u + r 2 r 3 (Au - 3y) 

+2x(y - u))d 3 + r 3 (rf - 2r 2 r\ + (r 2 + r 2 + 2u - x - y)n + r 2 (y - 2u))d 4 ), (B47) 
A\ = 4L2 "!f ' {-nr 3 T b dx + T b (n(4n - r 3 ) + 4(u - y))d 2 + T c (n(4n - r 3 ) + 4(u - y))d 3 - rir 3 T c d 4 ), (B48) 

A 2 = 4£ ^f 7/2 (-( ri - r 2 )r 3 T b rf 1 + T b (2x - 2r 2 - r 3 )d 2 + T c (2a; - 2r 2 - r 3 )d 3 - (n - r 2 )r 3 T c d 4 ), (B49) 

A 2 = 8L2 "!f ^ (T t (g - l)d 2 + T c (x - l)d 3 ), (B50) 

A\ = (r x r 3 d x + (n(4n - r 3 ) + 2(u - y))(d 2 - d 3 ) + nr 3 d 4 ), (B51) 

A 2 = -W /2 (ri-r 2 )^i+rf2-rf 3 + rf4), (B52) 
A 2 = 2 J L 2 m z 7/2 r 1 VF 3 "(di +d 4 ), (B53) 
A 2 = 2L 2 m z 7 / 2 (r 1 -r 2 )V^(d 1 +d 4 ), (B54) 

or 7/2 

A 2 = ,4 2 = - 2 ^_f (T b d 2 + T c d 3 ), (B55) 

^ii - ^12 = 0, (B56) 
m 9/2 N 

A i> = z ° ( r3Tb (rl _ r2ri + u + (_ 2 x - 2y + l)?/)di + T h {r 3 rl + (2ar(a; + y - 1) - r 2 r 3 )n + 2r 2 xy 

+r 3 (y + u - 2y 2 - 2xy))d 2 + T c ((r 2 r 3 + 2x(x + y- l))n - r 3 r\ + 2r 2 xy - r 3 (y + u - 2y 2 - 2xy))d 3 
+r 3 T c ((y - l)r\ + r 2 r x - u + (-r 2 + r\ + 2u)y)d 4 ), (B57) 

A% = -f ° (r 3 (y - r\)d x + (n(r 2 - r 2 n + 2u + g - 1) + (2r 2 - r 3 )y)ri 2 

+ (^r 2 - (r 2 + x - l)n + 2r 2 w - (2r 2 + r 3 )t/)d 3 + r 3 (r 2 - r 2 n + u)d 4 ), (B58) 
m z ^N 
4Li x /r 3 ~ 

+ (-2r 2 r 2 - 2yr 3 + r 3 + 2(r 2 - n)a;)d 3 - r 3 (r 2 - (n - r 2 ) 2 )d 4 ), (B59) 
Af - 9 /2 ^° (r 3 di + (2r 2 - r 3 + 4 n (x + y - l))d 2 + (2n - r 3 + 4r 2 (x + y - l))d 3 + r 3 (2a; + 2y - l)d 4 jiB60) 

Af = "! Z = ° {-rsTtyd! + r b (n(r 2 - r 2 n + u + a; - 1) + r 3 y)d 2 + T c {-r 2 r\ 



+ (r 2 + x + y - l)n - r 2 w + (r 2 + r 3 )j/)d 3 - r 3 T c (r 2 - r 2 n + u)d A ), (B61) 
mz 9 / 2 iVo, 

m z 9 / 2 7V , 



Aq = mZ ' '. (T b (2n(ri - r 2 ) + 2w + 2y - l)di + T h (2a; + 2y - l)d 2 + (T c - 2T c2 /)d 3 - 4 ri r 2 T c d 4 ), (B62) 

^1 



' (r 3 (-T 6 di + T 6 d 2 + T c d 3 ) + (2a; + 4y - l)(2T 6 nd 2 - 2T c r 2 d 3 - T c r 3 d 4 )), (B63) 



At = -^^y^(x + 2y-l)(T b d 1 +T c d i ), (B64) 



-4; 

A 
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A 4 ' 
A ii 



A 4 - 



m z 9 / 2 N 
m z 9 / 2 N 



2Li 
2Li 



(2nd 2 - (ri - r 2 )d 3 ), 
(nrf 2 + r 2 d 3 ), 
di, 
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(B65) 
(B66) 
(B67) 
(B68) 
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